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Abstract 

Flattened liquid drops poured on a very hot surface evaporate quite slowly and 
float on a film of their own vapour. In the cavities of a surface, an unusual type 
of vibrational motions occurs. Large vibrations take place and different forms of 
dynamic drops are possible. They form elliptic patterns with two lobes or hypotro- 
choid patterns with three lobes or more. The lobes are turning relatively to the hot 
surface. We present a model of vibrating motions of the drops. Frequencies of the 
vibrations are calculated regarding the number of lobes. The computations agree 
with experimental forms obtained in [1] by Holter and Glasscock. 
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1 INTRODUCTION 



Water splashed on a moderately hot metallic flat plate spreads out, comes to 
the boil and then quickly evaporates. It is not the same when the metal is 
very hot: water remains cool, breaks up into many drops that roll, bound and 
are thrown on all sides. Such drops are mentioned in the literature as being in 
" spheroidal state" . Most of the people have noticed these phenomena in their 
childhood when they stared at rolling drops on a very hot oven. 
Leidenfrost first experimented the process in 1756: a little water poured on a 
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red hot spoon, does not damp the spoon and takes the same shape as mer- 
cury. A strong motion of vapour lying between liquid and metal supports liquid 
masses and causes fast vibrating motion in the liquid bulk. This is the film 
boiling or Leidenfrost phenomenon [2]. It is well known since before the time 
of Faraday and was the subject of numerous studies during the nineteenth 
century [3]. 

One of the interests in the Leidenfrost phenomenon rises from an explosive 
behaviour associated with non-equilibrium of pressures. This affects all the 
industrial sectors where hot temperatures are used (for example in metal 
industry: quenching of metals [4], etc). This generates the unsteadiness of 
distillation plants in the petroleum industry and many accidents in nuclear 
engineering [5]. 

Indeed, the first surprising effect was previously described. The liquid drops 
are floating on a vapour film when they are placed on hot plates. This state is 
usually explained as the fact that when the temperature of the wall is greater 
than a value depending both on the fluid and the state of the surface, the 
exchanges of heat are small corresponding to the formation of a thin film of 
vapour isolating drops from the wall. In such conditions, liquid drops maintain 
during a time of the order of minutes and usual descriptions present the film 
of vapour keeping liquid drops well below boiling. Observations of the drops 
give phenomena by translational and rotational motions. Photographic and 
stroboscopic tools have been used to have descriptions of the experiments, 
but the effects can be seen by naked eye. These vibrations lead to beautiful 
patterns that appear as solid geometric figures. Such a phenomenon is quali- 
tatively very well described in the paper by Holter and Glasscock [1]. 
The different models do not give an interpretation of vibrating motions in 
the liquid bulk of drops on hot plates. It seems that classical treatments of 
vibrations (like in Lamb [6]) do not cover the present situation. Contrary to 
first impression, we determine that the present situation complicated by the 
gravity does not take into account heat gradients, vapour flows or more subtle 
effects. 

For the analytic approach, we consider a horizontal surface with a slight cavity 
invariant by rotation around its vertical symmetry axis. The drop is assumed 
not quite big enough to slip out of the cavity. Experiments describe such drops 
as flattened spheroids in the vertical plane. If the mean curvature radius of the 
surface is large compared to the size of the flattened drops, we assume that 
the motion of the liquid is approximately a motion in the horizontal plane. We 
consider plane motions of an incompressible liquid submitted to a convenient 
potential due to the effects of curvature of the surface. We investigate only 
rotations of drops and not random motions occurring by translation if the 
surface is really flat plate. 

Due to the film of vapour, we assume that motions are frictionless and the 
very small viscosity effect is neglected. The motions are not small vibrations. 
The form of drops in vibration are hypotrochoids [7] with different numbers 
of sides or lobes. The case of only two lobes corresponds to an elliptic form 
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of drops. The results of computation favourably compare to experiments. It 
is possible to evaluate the angular velocity of vibrations. It depends on the 
acceleration of gravity, the curvature of the surface and the number of lobes. 
Times of rotation for these spheroidal drops are computed. They seem to be 
in accord with simple observations [8]. 



2 A CLASS OF PLANE MOTIONS FOR INCOMPRESSIBLE 
PERFECT FLUIDS 

2. 1 Motions of a fluid 

The motion of a continuous flow can be represented by a surjective differential 
mapping 

C^x = HC), (i) 

where ( = (£, £) belongs to W, an open set in the time-space occupied by the 
fluid between time t\ and time ti- The position in the reference space T> is 
denoted by x\ its position at time t in V t is denoted by £ [9] . We assume that 
distinct points of the continuous fluid remain distinct throughout the entire 
motion. At t fixed, transformation (1) possesses an inverse. 

In order to describe a particular case of plane motions of a fluid let us introduce 
relatively to fixed Cartesian plane systems of coordinates, the inverse mapping 

Z^z = <p(t,Z). (2) 

The complex number Z = X + %Y is such that x — (X, Y) G 1Z 2 represents 
Lagrangian coordinates in the reference space V Q . The complex number z = 
x + iy with £ = (x, y) G TZ 2 represents Eulerian coordinates in the Euclidean 
plane T> t . Let us notice that (X, Y) is not necessarily a spatial position of a 
particle. Then, (X, Y) are material coordinates. 

In polar form Z = re l6 . Here we call (r,6) the reference position. Now, we 
write (r, 6) G V Q and in polar form, transformation (2) is written 

(r,6) eV ^ z = Vt (r,6) G V t . (3) 

Let us consider a particular transformation (2) of the form 

z = Z + J{Z)e^, (4) 

where a; is a real constant and the function f — p + iq is an analytic function 
of Z with p is the real part and q is the imaginary part of /. 
Components of the velocity V are (x,y) and z = iuf(Z) e' tuJt . 
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Components of the acceleration are (x, y) and z = ~uj 2 f(Z) e luJt . 

Let us denote F(Z) a primitive function of f(Z). At t fixed, relation (4) yields 



or 

and 
with 



e- tuJt f(Z)dz = e" tuJt f(Z)dZ + (p + iq)(dp - idq) 



e~ iuJt f(Z)dz = e~ iuJt dF(Z) + ]- d{p 2 + q 2 ) + i{qdp - pdq) 



TZ[e- luJt f(Z)dz} = d$ 

1 



$ = n[e-^f(Z)} + ^f(Z)f(Z) . (5) 
(71 denotes the real part of a complex number). 

Consequences of transformation (4) are: 

(a) xdy — ydx — i(xdx + ydy) = —i zdz = —u e~' tuJt f(Z)dz 

and consequently xdy — ydx — — cu d$ . 
Then 

div V = (6) 

and — oo $ is the stream function. 

(b) 'xdx + ydy + i(xdy — ydx) = 'zdz = —to 2 e~ wt f(Z)dz 

and consequently 'xdx + ydy = —u 2 d& . 
Then, 

T = -uo 2 grad $ . (7) 
Equation (7) shows that the acceleration V is the derivative of a potential. 

2.2 Pressure of these plane motions 



In the case of plane motion of an incompressible perfect fluid, the equation of 
motion yields 

T = -- grad p - grad W, (8) 
P 

where p is the constant volumic mass, p is the pressure and W is the extraneous 

force potential. 

From equation (7) we get 



p 
P 

Relation (9) yields the value of the pressure field as a function of $ and W. 



L = iU 2 $- W. (9) 
P 
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3 PARTICULAR CLASS OF PLANE MOTIONS 



For a reference set V Q such that 

V = [z = re ie with (re[O,r o ],0e[O,27r])}. 



(10) 



Let us consider two particular cases of motions defined by a function f such 
that: 

(a) f(Z) = 4^, neN, n > 1, aG^« 



(6) /(Z) = AZ, Ae]0,l[ 



5.1 Case (aj 
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is the representation of the motion. 

In parametric representation, relation (11) yields: 



(11) 



x = r cos 9 H cosiut — n9) 

y = r sin 6 H — — y sin(cjt — n#) 



(12) 



with (r, 9) E V . 



The trajectory of a particle with reference position (r,9) is a circle (C rj e). In 
the system coordinates, the centre of (C r ,e) is ( r cos 0, r sin 9) and the radius 
is R = r"/a n_1 . The particle associated (r,9) moves on the circle (C r ,e) with 
the angular velocity u>. 

At time t, particles whose reference positions in T> are on the circle with 
center O and radius r move in the physical space T> t on a hypotrochoid curve 
<£t(r). Curves <£t(r) are obtained as circular disks of radius p = r/n rolling 
internally inside a fixed circle of radius r + p. 

Two different points of V Q correspond to two different points of V t and the 
hypotrochoid has no double point. Due to 



det 



dx dx 

oy_ oy_ 

dr 89 



,2n-2 



r(l — n 



7 2n-2 



n-1 



the condition for the hypotrochoid has no double point is ( - I < — and 

\a) n 

the set V of material variables verify 



n— 1 | 
< - 

n 



(13) 
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Consequently, <£j(r ) is the free boundary of the fluid and the length a verifies 

a> r Q n -j/n . 

One verifies that 

<p t (r,8) = <p (r,e-nt)e 

with 



int 



n+1 

The hypotrochoids <£t(r) where r G [0,r o ] and specially the free boundary 
(£i(r ) of the fluid, are turning relatively to the hot plate with the angular 

velocity U = 



n + 1 
3.2 Case (b) 

Z eV ^z = Z + \Ze iuJt G V t (14) 
is the representation of the motion. 
In parametric representation, relation (14) yields: 

x = r cos 9 + \r cosiujt — 9) 

, ! (15) 

y = r sin 9 + Xr sin(ujt — 9) 

with (r, 9) G V . Contrary to case (a), no limitation is imposed on r Q . 
At time t, the particles whose reference positions in T> are on the circle cen- 
tered on O and radius r move in the physical space on an elliptic curve &t(r). 
Ellipses &t(r) have semi-axes of length r(l + A) and r(l — A). They are turning 
relatively to the hot plate with the angular velocity f2 = uo/2. 



4 MOTIONS OF A PERFECT FLUID IN THE BOTTOM OF A 
LARGE CAVITY. 



With a convenient Euclidean frame Oxyz, the equation of the hot surface is 

z = f(r) with r 2 = x 2 + y 2 . (16) 

Function / is a C 2 -function and /(0) = f'(0) = 0. The vertical direction is 
Oz. Equation (16) yields 

f"(0) 

z = —^-r 2 + r 3 e(r) with lim r ^ e(r) = . (17) 

For a surface whose meridian curve is of small curvature, we use the equivalent 
approximation 
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The potential of gravity forces is 



W = k(x 2 + y 2 ) with k — g 



f(0) 



(19) 



Relation (18) is used in the following computation: the mean curvature of the 
surface at (0, 0, 0) is it! = and 



k 



_9_ 

2R 



Let us consider the two cases (a) and (6): 
In case (a) 

Z n 



7n+l 



f(Z) 



,n- 1 



F(Z) 



(n + 1) a™- 1 



We deduce the stream function, 



„2n 



,n+l 



2a 2 ™- 



For the potential (19) we deduce from relation (12) the expression 



W = k 

Equation (10) yields 



„2n 



„n+l 



,2n-2 



+ 2-^- cos (cut - (n + 1)6>) 
a 



(20) 



„2n 



„n+l 



-2a 2n ~ 2 (n + l)a 

^2n ».n+l 



— cos(c^- (ra + l)0)] 



Mr 2 + + 2— T cosM - (n + 1)5)] + C7 te . 



In the case u> 2 = 2(n+l)/c and for a pressure p on the free boundary associated 
with r = r Q , we obtain 



P Po 



UJ 



p p 2(n+l) 
In the case (6) 



irl ~ r 2 ) - 



n 



UJ 



yi 2??. ,v! 2t?< 

2 'o ' 



2(n + l) a 



2n-2 



(21) 



and 



A A 2 A 

F{Z) = -Z 2 , $ = — r 2 + -r 2 cos(c^-2#) 
2 2 2 



W = A;r 2 [l + 2A cos(ut - 29) + A 2 ] . 
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Equation (12) yields 
A 2 A 

V - = r V [— + - cos(ut - 26)} - kr 2 [1 + 2A cos(ut - 26) + A 2 ] + C te . 
In the case u 2 = 4k and for a pressure p a on the free boundary, we obtain 

2 

- = - + T (l"A)(r -r ). 

In the two cases, equation (20) yields u 2 = (n + l)g/R and fixes the angular 
rotation of the drop 

n 2 = - — g —— . (22) 

Consequently, following the different values of n, different modes are possible. 
The period of the rotation of the vibrating drop is given by 



T = 2K^ n + 1)R . (23) 

In fact, after an interval of time T n+1 =T/(n+ 1), the drop is identical to its 
previous position and v n +i = l/T n+ i is the apparent frequency of changes. 



1 (n + l)g 



5 COMPUTATION AND GRAPHS OF THE MOTION 



Experiments were performed with a very simple apparatus in [1]. No data were 
performed, but it seems possible to investigate the above frequencies. Some 
apparent frequencies of the drop vibrations are calculated in Table 1. 
We have drawn different cases of vibrating drops from two lobes (basic mode) 
to eight lobes. Drops are turning around their center of symmetry and we 
have drawn different positions of a drop with two and three lobes (Fig. 1). 
In the case of three-lobed modes, the trajectory of a particle is represented 
as a circle for a particular drop. The graphs have the same forms than the 
ones given in [1]. Only the case of two lobes present some discrepancy with 
the interpretation given in [1]. In [1], the middle of the drop may be, in some 
cases, squeezed by comparison with the shape given here. Depending on the 
value of a and n, the drop can be like vibrating polygons of any type (Fig. 2). 
The two orientations are similar by computation; we notice also it is possible 
that the drop reverses the direction of rotation of the motion. 
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30 


1.29 


1.58 


1.82 


2.04 


2.23 


40 


1.11 


1.37 


1.58 


1.76 


1.93 


50 


1.00 


1.22 


1.41 


1.58 


1.73 



Table 1 

Some apparent frequencies 7y n+ i (in Hertz) of drops as function of the mean 
curvature radius of the surface (in cm) and the number n+1 of lobes of the drop 
(g=981 cm/s 2 ). 
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Two- lobed mode 




Three- lobed mode 



Fig. 1. In the case of two-lobed mode or three- lobed mode, the drops are turning 
around their centre of symmetry. The trajectory of a particle is a circle represented 
in the case of a three-lobed mode. 
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Fig. 2. Two possible forms of drop are presented in cases of three-lobed to eight-lobed 
modes. In many cases the vibrating drops are approximately polygons which may 
have any of three to n sides. 



